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Abstract. 

We review recent studies about the resonance spectrum of quantum scattering 
systems, in the semiclassical limit and assuming chaotic classical dynamics. 
Stationary quantum properties are related to fractal structures in the classical 
phase space. We focus attention on a particular class of problems that are chaotic 
maps in the torus with holes. Among the topics considered are the fractal Weyl 
law, the formation of a spectral gap and the morphology of eigenstates. We also 
discuss the situation when the holes are only partially transparent and the use of 
random matrices for a statistical description. 



1. Introduction 

A closed quantum system has a discrete spectrum of energy levels, to which are 
associated bound states. Open systems, on the other hand, have a resonance spectrum, 
consisting of a discrete set of complex poles of the Green's function or of the scattering 
matrix (besides having a real continuous spectrum of scattering states) . If we imagine 
a closed system being slowly opened, its bound states will acquire small (negative) 
imaginary parts and become resonances. In this sense, the resonance spectrum can also 
be seen as eigenvalues of a non-hermitian Hamiltonian. If £ = E — iT/2 is a complex 
energy value, the usual time evolution factor \ e ~ l£t / h \ 2 = e -rt/h j g n0 ^ un i mo dular 
but decreases exponentially. The quantity T/K is thus called the decay rate of the 
state, or equivalently h/T is called its lifetime. Insofar as any real system is always in 
contact with its environment, most energy levels are actually resonances with a finite 
decay rate. 

We are interested in systems for which the classical dynamics is chaotic. 
Scattering of waves in chaotic systems has been experimentally realized in a great 
variety of physical systems. Out of the vast literature available, we select only a few 
examples: electron transport in semiconductor quantum dots pQ [21 El H] and graphene 
quantum dots [SJ [7] ; microwaves in normal metallic [SI [SJ [TO] or superconducting 
cavities [IT] ; nuclear reactions [12j [13] ; acoustic waves [14] ; microlasers [15] [SI [17] . 

In the semiclassical limit H — > the wavelength is much smaller than any classical 
length scale and the nature of the classical (or ray) dynamics becomes important. 
Three main questions then stand out in connection with the resonance spectrum of a 
chaotic scattering system. First, how does the number of resonances with given decay 
rate grow with ft? Second, is there a lower bound for the decay rates? Third, what do 
the resonance wavefunctions look like (inside the scattering region)? As we will see, 
all these questions are related to fractal properties of the classical dynamics. They 
are all still open to some extent, although much progress has indeed been made. 
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The systems we have in mind have 'classical' openings. By this we mean the 
following. For a fixed value of the energy, a finite 'energy shell' can be constructed by 
imagining a big enough box enclosing the scattering region in configuration space and 
all possible momenta with fixed magnitude. In this energy shell, the points leading to 
immediate escape occupy a finite volume. In this situation the resonances are always 
strongly overlapping (sometimes called the Ericson regime). 

A very popular class of models with which one can theoretically study chaotic 
scattering are so called torus maps, some of which are well known like the baker map 
or cat maps. These toy models have the maximal simplicity still allowing for chaotic 
dynamics. Also, they can be quantized in order to study wave properties. In the last 
few years, the resonance spectrum of open quantum maps has been widely studied and 
have shed considerable light on the more general problem. Our purpose is to review 
these developments. 

Other reviews concerning areas related to the ones discussed here have appeared 
recently and can offer complementary insights. Open systems and non-hermitian 
Hamiltonians were reviewed by Rotter in [18] and, from a very different point of view, 
by de la Madrid and Gadella in [19] (see also [20]); open billiards were reviewed by 
Dettmann in 2T| and also by Altmann, Portela and Tel in [22]; the random matrix 
theory (RMT) approach to resonances in quantum systems was reviewed by Fyodorov 
and Savin in [53] and by Mitchell, Richter and Weidenmuller in [21] . Finally, spectral 
results for chaotic open maps were reviewed by Nonnenmacher in [25) . but from a 
mathematically more rigorous perspective than the one adopted here. 

In Section 2 we introduce open chaotic maps. We give a few examples and 
discuss their long time behaviour. We introduce fractal trapped sets and conditionally 
invariant measures. In Section 3 we discuss open quantum maps and their general 
properties. Section 4 is dedicated to the fractal Weyl law that has been conjectured to 
hold for the resonance spectrum in chaotic scattering. Section 5 discusses the possible 
existence of a spectral gap. The phase space morphology of resonance wavefunctions 
is considered in Section 6. In Section 7 we focus on systems where escape is not 
ballistic, like e.g. microlasers. The RMT approach to quantum scattering is discussed 
in Section 8. Finally, we close in Section 9 with some open problems. 

Let us mention that we do not consider the closely related, and somewhat 
complementary, transport formulation of scattering [26], which deals with S matrices, 
transmission eigenvalues, counting statistics, etc. The semiclassical and the RMT 
approaches to quantum chaotic transport were reviewed, for instance, in [27j and [28], 
respectively. 

2. Classical maps 

When studying the integrability of the dynamics of the solar system, Poincare 
introduced the idea of 'cutting through' phase space with a hyperplane, and recording 
the intersections of the system with this plane. This construction, now known as a 
Poincare section, produces a discrete-time dynamics in which the trajectory becomes 
a sequence of points. The time interval between intersections is lost, but this is usually 
considered a small price to pay for the dimensional reduction achieved. 

We consider maps in themselves, without any reference to a continuous-time 
system from which they could be derived. We shall treat only maps that are defined 
on a two-dimensional torus, which can be represented simply as a square with opposite 
sides identified. Let us denote by {q,p) orthogonal coordinates in this space which are 
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the sides of the square. We have in mind conservative chaotic dynamics (hyperbolic 
and mixing). 

Let M. denote the map. If (q',p') = M(q,p), we say that the point (q,p) evolves 
into the point (q',p'). An infinite sequence of points is called an orbit. If an orbit 
consists in the infinite repetition of a finite number of distinct points, it is called 
periodic. The number of distinct points is then its period. Its is known that chaotic 
systems have infinitely many periodic orbits, and that they form a dense set in phase 
space, i.e. there is a periodic orbit arbitrarily close to any point. 

Many books discuss chaotic maps. We just mention for example 29, 30], where 
all basic definitions, examples and further discussion can be found. In the following, 
we briefly present a few paradigmatic systems and introduce only the notions that will 
be needed when we consider the resonances of open quantum maps. 

Perhaps the simplest example is the baker map, given by 

, , f(2g,p/2), if q < 1/2, 

( ^l( 2g -l l(p+ l)/2), if, > 1/2. (1) 

Its well-known dynamics consists in a stretching by a factor of 2 along the q axis, 
together with a contraction by the same amount in the p axis, plus a 'cut and paste' 
operation involving the line q — 1/2. Clearly, area is preserved. Also, the map is 
hyperbolic, with stable and unstable directions at every point being parallel to the 
coordinate axes. The Lyapounov exponent is In 2. Notice that the baker map may be 
written as 

(?,p).->(2g-[2?],(p+[2?])/2), (2) 

where [2q] denotes the integer part of 2q. 

Many generalizations of the baker map exist. One of them makes use of D basic 
regions, stretching and compressing by a factor of D: 

(q, P )^(Dq-[Dq],(p+[Dq])/D). (3) 

This is sometimes called the D-nary baker map. The usual baker map corresponds to 
D = 2. This map has a Markov partition with D cells, and each cell is stretched and 
compressed by a factor D at each iteration. Therefore, the Lyapounov exponent in 
this case is InD. 

The so-called cat maps provide another paradigmatic family of chaotic linear 
maps. The following is an example: 

{q,p) ^ {2q+p,3q + 2p)modl. (4) 

In this case the Lyapounov exponent is In (2 + y/S). The dynamics of cat maps has 
no discontinuity. The stable and unstable directions are still the same at every point, 
but they are no longer orthogonal. 

Our final example will be the standard map. It is defined as 

(q,p) i ^ (q + p + Ksm(2nq),p + Ksm(2nq)). (5) 

It is more generic than the previous examples because it is non-linear. Its dynamics is 
known to be predominantly chaotic if K > 7. The stable and unstable directions, as 
well as the local stretching factor, are no longer constant in phase space. It is known 
that the Lyapounov exponent (the average stretching factor) is given approximately 
byln(K/2). 

Instead of propagating points, it is often useful to consider propagation of 
probability measures. Given a probability measure fx, its evolution is another 
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probability measure // such that, for any subset A of phase space, n'(Ai(A)) — fi(A). 
A measure /i is called invariant if fi(A4(A)) — (i(A) for every A. If we allow singular 
invariant measures, such as linear combinations of Dirac deltas, then clearly there are 
infinitely many of them (for example, to any periodic orbit we can associate one). 

If a probability measure fj, has a probability density /, such that fJ,(A) = 
f A f(q,p)dqdp 1 it is called absolutely continuous. Probability densities can also be 
evolved in time: if (q',p') = -M.(q,p), then the function g(q,p) such that g{q',p') = 
f(q,p) is called the evolution of / (remember that the map is conservative, so there is 
no Jacobian). In other words, each point simply carries along the value of the function 
associated with it. The linear operator V that implements this evolution, g = Vf, is 
called the Perron- Frobenius operator associated with the map A4. 

All examples of chaotic maps discussed previously are ergodic systems and hence 
have only one invariant probability density. This is called the natural density or 
equilibrium density, denoted f e (for conservative maps this is simply the constant 
function). They also have the exponential mixing property, which implies that any 
smooth initial density function will converge to f e exponentially fast in time (when 
defined on an appropriate function space, V has a single eigenvalue equal to 1 and 
all other eigenvalues have strictly smaller modulus). This convergence of course takes 
place in a weak sense. 

We now wish to consider open maps as models of chaotic scattering. Given a 
map, an open version is obtained by defining a 'hole' in phase space, which can in 
principle be any region of finite area (or union of regions) . Once a point falls into the 
hole, it stops being propagated (i.e. it 'escapes'). This means that the dynamics is no 
longer conservative. In fact, for chaotic systems almost all initial conditions eventually 
escape because of ergodicity. We denote by A4 the open map, i.e. the map that acts 
by first removing points in the hole and then evolving the remaining ones according 
to M. 

The trapped set of the dynamics is the set of initial conditions which do not 
originate from the hole and never hit the hole. In other words, they remain in the 
system for infinite time, whether propagated forwards or backwards. This invariant 
set will be denoted Kq. It is quite common in the quantum chaos literature to refer to 
Ko as the system's 'rcpcllcr'. This has been criticized on the grounds that 'repeller' 
should be used for sets that are unstable in all directions. When a stable manifold 
exists, [22] and [31] suggest calling K the chaotic saddle or simply the 'saddle'. We 
shall follow this suggestion. 

The stable manifold of the saddle, the set of points which converge to it and 
therefore never escape in the future, is also known as the forward-trapped set and 
denoted K + . On the other hand, the unstable manifold contains the initial conditions 
which never escape in the past and is also called the backward-trapped set K_. 
These two sets are fractals and have similar structures. Globally, they arc very 
convoluted. Locally, they are continuous in one direction and fractal in the orthogonal 
one. Therefore, they both have zero measure in phase space. 

The sets K- and K + can be described in terms of images and pre-images of the 
hole. Let h m — A4 m h be the mth image of the hole under the dynamics. The set 
h° = h is the hole itself.J^hese sets can have finite intersections, so we define another 
family of sets as lZ m — A4 m h. The set lZ m contains the points that fall into the hole in 
exactly m steps when propagated backwards. Analogously, we can start with 1Z = h 
and for m > define Hr m as the set of points which fall into the hole in exactly m 
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Figure 1. (color online) Important sets for the open ternary baker map. The 
region 1/3 < q < 2/3 corresponds to the hole. The set of points that fall in 
the hole in exactly one step (IZ^ 1 ) have q in the region (1/9, 2/9) |J(7/9, 8/9). 
The darkest lines in the left panel are an approximation to the forward trapped 
set, K+. In the middle panel we see the sets TZ 1 as 1/3 < p < 2/3, 1Z 2 as 
p 6 (1/9, 2/9) U(7/9, 8/9), etc. The darkest lines in this panel approximate K-. 
The right panel shows a finite-time approximation to the invariant saddle Ko, 
which is the product of two Cantor sets. Figure reprinted with permission from 
[32]. 




Figure 2. (color online) Similar to the left panel of Figure 1, but for the kicked 
rotator. Different colors (shades) correspond to points that escape in one, two, 
three and four steps. The white region is a finite-time approximation to K+. 
Figure reprinted with permission from 1701 . Copyright (2004) by the American 
Physical Society. 

steps when propagated forwards. These sets obey AAlZ m D lZ m+1 , and 

MK m =K m+1 . (6) 

Let us imagine that we start with the torus completely filled with initial 
conditions, and they are all propagated forwards. First, the points at the hole escape. 
Then, the points that were initially at TlT 1 escape. Next, the points initially at 1Z~ 2 , 
and so on. Clearly, since it contains the points that never escape, the forward trapped- 
set K + must be the complement of the union of all lZ~ m . On the other hand, the 
backward-trapped set K_ is the complement of the union of all lZ m . Notice that the 
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sets lZ~ m intersect the unstable manifold K_ and in fact provide a useful partition of 
it (every point of K_ is either in one of the lZ~ m or belongs to K ). 

The simplest example of open map with a nontrivial saddle is the ternary baker 
map with the middle region identified with the hole. This is discussed in detail, for 
example, in [32]. It is easy to see that all sets 1Z ~ m consist of vertical strips. The 
horizontal section of TZ~ X , for example, is (1/9,2/9) U (7/9,8/9). The set K~ m is 
always the union of 2 m strips of width 3 _m . As a matter of fact, the forward trapped 
set is the product of the interval [0, 1] in the vertical direction and a Cantor set in 
the horizontal direction. Analogously, the backward trapped set is a Cantor set in 
the vertical direction. Finally, the saddle Kq is the product of two Cantor sets. This 
example can be seen in Figure 1, while the analogous construction for the kicked 
rotator is shown in Figure 2. 

Open systems can have invariant measures, but they must be supported on the 
saddle and thus are very singular. Let /i* denote the i-step propagation of measure /z. 
A measure satisfying /i* = e~ Tt fi is called an eigenmeasure of the system. Invariant 
measures correspond to T = 0. 

Incidentally, open systems provide a glimpse at how truly vast is the set of 
invariant measures for closed ones. A closed system can be opened by choosing any 
region of phase space as the hole. In turn, this may generate a chaotic saddle, which 
can carry many invariant (singular) measures. These are also invariant measures of 
the closed system, and we thus have many such measures for many possible choices of 
hole in phase space. 

Eigenmeasures with < T < oo are called conditionally invariant measures 
(in the sense that they are invariant provided we renormalize them at every step). 
Conditionally invariant measures for maps were reviewed in |33| . They are necessarily 
supported on K-. Since this set has a fractal 'cross section', they cannot be absolutely 
continuous. If we take smooth initial measures and propagate them forwards, they all 
still converge (when renormalized) to the same measure fi e , called the natural measure 
or equilibrium measure. It is an eigenmeasure satisfying 

= e" 7 V (7) 

The quantity 7 appearing in ([7]) is called the system's decay rate. Decay rates for 
chaotic systems with holes have been reviewed in |22j . In the limit of vanishingly small 
holes, 7 is given simply by the size of the hole. For finite holes, however, it is in general 
a complicated function of its size, shape and position (see the many references in [22 ). 
For example, the decay rate can be reduced if the hole contains a short-periodic orbit, 
because this implies that it overlaps considerably with its images, reducing the amount 
of points that escape after one step (this was observed for example in [34, 35, 36il37]). 

A formula for 7 is available in terms of the equilibrium measure |38j : 

7 = -ln(l- Me (/i)) ) (8) 

where fi e (h) is the /i e -measure of the hole h. This can be understood as follows: 
for long times any initial probability measure converges to )i e , and fi e {h) gives the 
proportion of points which escape at one time step, which must be equal to 1 — e~ 7 . 
It is also possible to compute 7 using a sum over periodic orbits contained in the hole 

[SSIEQ]. 

One important characteristic of a fractal set is its dimension. Because K- is a 
simple line in one direction, its total dimension is given by 

d(K-) = l + d(K_), (9) 



Resonances in open quantum maps 



7 



where d is called the partial fractal dimension. If the dynamics has time-reversal 
symmetry, then K + has exactly the same dimension as K_. Their intersection, K a , 
is locally equal to the product of their fractal parts and thus has 

d(K ) = 2d(K.) (10) 

as its dimension (for systems without time-reversal symmetry, the dimension of Kq is 
the sum of the partial dimensions of K— and K + ). 

We only consider fractal sets (K-, K + and Kq) characterized by a single fractal 
dimension (i.e. their Minkowski and Hausdorff dimensions coincide). For fractal 
measures, on the other hand, different notions of dimensions coexist. The most used 
ones are the Minkowski (or box-counting) dimension do, the information dimension d\ 
and the correlation dimension d 2 ■ These are in principle distinct but in most practical 
situations they happen to be quite close to each other. It is known that they can be 
generalized to a continuous function dp which is non-increasing (the so-called Renyi 
dimensions), so do > d± > d%. 

We are only interested in the natural measure /i e . Since it is supported on 
its Minkowski dimension is given by do(/j, e ) = d(K-). Since it is continuous in the 
direction of K-, we can write, for every /3, the corresponding dimension in terms of a 
partial component: 

^(/i e ) = l + ^(/i e ). (11) 

We will use d without a subscript for the fractal sets and with a subscript when we 
refer to dimensions of the natural measure. It is quite common, however, to talk about 
dp{K-) when one means dp{fjL e ), because of the uniqueness of the natural measure. It 
is well known that there exists an important equation (the Kantz-Grassberger relation 
|39j ) involving the decay rate, the Lyapounov exponent and the partial information 
dimension of /i e . This is 

dM = i~l- (12) 

Coming back to the open tribaker map, if we fill phase space uniformly with initial 
conditions, then exactly one third of points will hit the hole at each step, leading to 
7 = — ln(2/3) = In 3 — In 2. The partial dimension of K_ is the fractal dimension of 
the Cantor set, well known to be d(KJ) = In 2/ In 3. The natural measure is constant 
on K_ and, therefore, all its fractal dimensions are the same, dp(fi e ) — d(K^). This 
is not a general feature; it comes from the fact that the stretching rate is constant for 
this system. Notice that the Lyapounov exponent is A = In 3, so that the relation (|T2"j) 
is indeed verified. 

In Section 7 we shall discuss maps for which escape is not ballistic, i.e. when a 
point falls in the hole, it has a certain probability of escaping but can also continue 
being propagated. A probability density is then only partially attenuated when coming 
in contact with the hole, instead of being drastically cut. This is supposed to model, 
for instance, refractive escape of light rays from a dielectric sample. 

3. Quantum maps 

Maps have long been used as toy models in classical mechanics, because they allow for 
chaotic conservative dynamics even in a two-dimensional phase space. But they have 
also enjoyed popularity in the context of quantum mechanics, the reason being that 
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their numerical implementation is quite simple. When a map is defined on a compact 
phase space, for example, its quantum propagator is automatically finite dimensional. 

Unlike for Hamiltonian systems, there is no standard procedure to quantize a map. 
The general idea is that, given a conservative classical map M, one should associate 
with it a unitary operator U (the propagator) acting on some Hilbert space. Naturally, 
some kind of classical limit must be defined, in such a way that the dynamics of A4 is 
recovered. 

Since we are considering only the torus as our phase space, we must enforce 
periodicity both in q and in p. This implies that both directions become discretized, 
i.e. there will be a finite number N of 'position' eigenstates \q n ) and 'momentum' 
eigenstates \p m )- Quantization of classical periodicity can be achieved with a pair of 
arbitrary phases, so that 

(q n + = (<Z„|V)e 2 ™ X «, ( Pn + = (Pn\lp)e 2 ™*<>. (13) 

For given values of n and m, the corresponding coordinates on the torus are 

n + X P m + x q n .s 

The case Xg = Xp = corresponds to periodic boundary conditions, while Xq — Xp = 
1/2 correspond to antiperiodic boundary conditions. Position and momentum must 
be Fourier related, 

(Pm|9«) = (F^ Xp ) nm = le- 2 "(» + ^^«)^. (15) 



Here both n and m range from to N — 1. The dimension N takes the role of inverse 
Planck's constant: 

N =^-T- (16) 
2irn 

Quantization of the dynamics consists in the construction of a quantum 
propagator U, which is a iV-dimensional unitary matrix responsible for taking one 
wavefunction to another: 

N-l 

Tp'iQn) = ]T U nk i>(q k ). (17) 

fc=0 

This operator can always be diagonalized, and its spectrum lies in the unit circle in 
the complex plane, namely all N eigenvalues are of the form e i6n with < 6 n < 2tt. 

The matrix U must somehow reproduce the classical dynamics in the limit 
N — > oo. To visualize this limit one usually resorts to tools such as coherent states and 
Husimi functions. In usual quantum mechanics, coherent states \q,p) are ground states 
of harmonic oscillators, i.e. Gaussian wavepackets with minimum uncertainty. When 
realized in the torus, they must incorporate periodicity (a recent review about coherent 
states can be found in [40]). Therefore, they are defined in position representation as 

oo 

{q n \q,p)=C e- 2 "^e-^('»-' +t » 2 e 2lW («»-«/ 2+fc ) ,, 1 (18) 
k=— oo 

where C is a normalization constant. 

The Husimi function %^{q,p) — \{ij)\q,p)\ of a quantum state \tp) is a real, 
non- negative function defined over all phase space. It can thus be normalized and 
interpreted as a probability density (however, it does not have the correct marginals, 
e.g. J T-L^(q,p)dp 7^ \ip(q)\ 2 ). The Husimi function of a coherent state, W\Q,p)(q,p), 
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is a periodized Gaussian in both directions, centered at the point q = Q,p = P, with 
widths ~ 1/VN. 

For large N, the action of the quantum propagator on coherent states must satisfy 
the semiclassical condition U\q,p) ~ \A4(q,p)}, i.e. it must take the point q,p to 
its image under the classical map. Actually, this evolution will always introduce a 
stretching of the wavepacket along the unstable direction, coming from the Lyapounov 
exponent. It is more correct to say that the action of the quantum propagator must 
correspond, in the semiclassical limit, to the action of the Perron-Frobenius operator: 
it takes an initial probability distribution centered at (q,p) (the Husimi function of 
the state \q,p)) to another probability distribution which is centered at Ai{q,p) (but 
is not the Husimi function of the state \M{q,p))). 

This semiclassical approximation cannot hold for long times due to interference 
effects. The similarity between the action of U* on the quantum side and V* on the 
classical side must degrade with time t. For chaotic systems, it degrades exponentially 
fast. Heuristically, it will break down at the time scale given by the time it takes a 
minimum wavepacket (of width ~ 1/ \/N) to stretch to the size of the system. This 
is called the Ehrenfest time and, since the stretching is regulated by the Lyapounov 
exponent according to e At , it is given by InN. 

Notice that other definitions of Ehrenfest time are possible. For example, the 
time it takes to stretch a state which is squeezed along the stable manifold into a state 
which is squeezed along the unstable one is jhiiV. This is equivalent to stretching 
the smallest length, 1/N, to the size of the system and has also been called the 
Ehrenfest time in the literature (indeed this is the time when some interference effects 
start to become important, see for example Til). Yet another notion of Ehrenfest 
time, related to escape through holes, will be presented below. The quantum-classical 
correspondence for open chaotic systems, with special emphasis on different Ehrenfest 
times, was reviewed in |42j . 

Let us now mention, without any technical details, some of the quantizations that 
have been devised for the classical maps discussed in the previous section. 

The first quantization of the baker map, introduced in [43] . postulated the 
following propagator: 



It is compatible with periodic boundary conditions and can only be defined on even- 
dimensional spaces, reflecting the nature of the underlying map. It can be verified 
that this operator has the baker map as its classical limit, but an inversion symmetry 
of that map does not hold exactly. Later |41) . another quantization with antiperiodic 
boundary conditions was suggested, that had exact symmetries. In this quantization 



The quantization of cat maps was considered in [44l [45j US] . Here the fact that 
the classical action is quadratic was used to define a semiclassical quantization. For 
the particular cat map Q we mentioned previously, this results in a matrix U whose 
elements are quite simple: 




(19) 



Fjj is replaced by F } 



1/2,1/2 



N 




(20) 



with n, m between and N — 1. 
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Finally, the standard map can also be quantized [J71 H8]. It can be realized 
as a kicked rotator, which leads to a quantization in terms of the product of free 
propagation and kicking. More convenient is the following closed formula in position 
representation: 



This actually corresponds to a version of the system which displays time-reversal 
symmetry. 

A few experimental realizations have been suggested for the quantum maps 
discussed above. In [49] an optical set-up of the baker map was proposed, while a 
possible realization using nuclear magnetic resonance was discussed in [50 . Already 
in [H] the possibility of realizing the quantum cat map using Fresnel diffraction by a 
periodic grating was suggested. Finally, the quantum kicked rotator has already been 
realized in the laboratory by placing cold atoms in a pulsed standing wave [5TJ [52] . 

It is natural to expect that in the semiclassical limit the Husimi functions 
of eigenstates of U should converge, in some weak sense, to invariant probability 
distributions. As we have seen, these are either related to f e or singular. It is known 
that, as N — > oo, almost all eigenstates will become equidistributed, i.e. their Husimi 
functions will converge to f e . In other words, as one diagonalizes U in larger and 
larger dimensions, the fraction of equidistributed states should grow to 100%. This 
property is called quantum ergoditicy. A stronger property is that of quantum unique 
ergodicity. This property holds when we can say that all quantum states become 
equidistributed without exception, which is a stronger statement than saying that 
their fraction becomes 100%. Quantum ergodicity and quantum unique ergodicity 
have been extensively investigated for cat maps and baker maps (see, for example, 
[56] 157] 158] [59]). A recent review appears in [60] . 

Numerically, it is easy to find for finite N some eigenstates with Husimi functions 
concentrated on periodic orbits. This phenomenon has been termed 'scarring'. For 
continuous-time systems, there is a large literature on this subject starting with |53j ; 
for maps see [54] [55] and also [60] . 

In order to open a quantum map, what we have to do is identify a sector of 
the Hilbert space with the hole. This is done by 'projecting out' of it, i.e. the open 
propagator is 



where U is the closed propagator and n is the projector on the complement of the hole. 
It is most common to use a hole that is a strip of width 8 parallel to the momentum 
axis. In that case n will be diagonal in position representation. The value of 6 sets a 
new length scale in the problem [42]. We can define the escape Ehrenfest time r e as 
the time it takes to stretch the fundamental size 1/N up to 5. This satisfies e AT " = NS. 
If we let K be the number of quantum states that fit in the hole, i.e. the dimension 
of the kernel of the projector II, then clearly r e = i In if. 

The eigenvalues of U will be located inside the unit disk in the complex plane. 
They can be written as 

z n = e^~ r "/ 2 , (23) 

so that r n can be interpreted as a decay rate (it corresponds to what was T/H in 
the more general Hamiltonian setting). The set of all N eigenvalues comprises the 




U = UIL, 
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Figure 3. (color online) Numerical verification of the fractal Weyl law. In the 
left panel [Figure reprinted with permission from [71] : Copyright (2008) by the 
American Physical Society], the logarithm of the number of states having decay 
rate smaller than some constant value is plotted against InN, for the kicked 
rotator with different holes. In the right panel [Figure reprinted with permission 
from |72[ : Copyright (2008) by the American Physical Society], the fraction of 
states having decay rate smaller than some constant value is plotted against N, 
for the cat map with different holes. In all cases, the FWL prediction (lines) is in 
excellent agreement. 



resonance spectrum of this 'scattering' problem. Resonances in the vicinity of the 
origin as seen as short-lived, while any finite value of r„ can be understood as being 
long-lived. 

4. Fractal Weyl law 

The Weyl law is a very basic result in the asymptotics of energy levels for closed 
systems. In a nutshell, it says that each stationary quantum state occupies a volume 
(2ttH) 1/ in phase space, where v is the number of degrees of freedom [61]. This means 
that, to leading order in h, the number of states with energy less than E is proportional 
to fr v . Alternatively, let fl(E) be the set of phase space points with energy E. It is 
obvious that this set has dimension do = 2v — 1. The density of states around E is 
then proportional to the size of Vt(E) and scales with h as fi - ( d °+ 1 )/ 2 . 

For open systems, a similar idea can be applied to the resonance spectrum. One 
can choose a small but fixed region of the complex plane, around a point with real 
part E and imaginary part —hT/2 (such that the decay rate is fixed), and ask how 
the number of resonances inside it grows as h —¥ oo. For continuous-time chaotic 
flows, it is conjectured that this number should behave like h~^ do+1 ^ 2 , where now do 
is the fractal dimension of the classical chaotic saddle at energy E. This is known as 
the fractal Weyl law (FLW). In the context of continuous time, it has been checked 
for scattering from three gaussian 'bumps' [53] [53], three rigid disks [53], four rigid 
spheres in 3D [55] and a modified Henon-Heiles potential [66] . The only rigorous result 
available to date is that the exact number of resonances is bounded from above by the 

FWL [62 Eg. 

A vague heuristic justification of the FWL can be presented as follows (a more 
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Figure 4. Numerical verification of the fractal Weyl law for the baker map. 
The quantity n(r, N)/N~ d / 2 is plotted against r, where n(r, N) is the number of 
eigenstates of the map with modulus smaller than r and d is the fractal dimension 
of Kq. The plot shows that the results for many values of N lie on top of each 
other, validating the scaling. From |75l . reprinted with permission. 



detailed account is given in 62 ). States with energy E — iT/2 have lifetime H/T. From 
basic principles, these states should be related to structures in classical phase space 
which remain inside the system for that long. By analogy with the usual Weyl law, 
their number is estimated as being proportional to the volume of a 'coarse-grained' 
version of the saddle, of thickness proportional to h. As h — > 0, this volume is related 
to the fractal dimension. 

The FWL takes on a particularly simple form for quantum maps, because of 
dimensional reduction, no energy-dependence of the dynamics and the fact that the 
total number of resonances is equal to N. In that case the conjecture is that the 
number of resonances with modulus larger than some r should grow like C(r)N d / 2 , 
where C(r) is some unknown function of r and d is the fractal dimension of the saddle 
Kq discussed in Section 2. Notice that 



where d(K_) is the dimension of the backward trapped set and d{K-) is the associated 
partial dimension. In [32j [69] the FWL was verified numerically for the open triadic 
baker map and rigorously proven for a non-canonical quantization of this map. It has 
also been verified for the kicked rotator [701 E] an d f° r the cat map [75] . Again, only 
upper bounds can be rigorously proven |73j . In Figures 3 and 4 we show numerical 
verification of the FWL for three types of open quantum maps. 

What is the function C(r)? The first point that needs to be settled is whether 
this function is universal or system-dependent. For the kicked rotator [70] it is quite 
close to the universal prediction from random matrix theory (see Section 8), but this 
is not true for the baker map [32]. In fact, in this latter case there seems to be a gap 
in the middle of the spectrum, for which no explanation has been offered so far. If the 
function C(r) is system-dependent, it remains to be determined what kind of classical 
information goes into it. 

One theoretical approach to the FWL, put forth in [70], is via the construction of 
short-lived states. We sketch its argument. Let K be the dimension of the kernel of 



d(K ) 



d(K-) - l = d(K^), 
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Figure 5. (color online) Illustration of how propagation up to Ehrenfest time 
can lead to almost perfect escape of a wave packet. The initial packet, whose 
center lies in 1Z~ S , is very stretched along the stable manifold of the dynamics. 
Upon propagation, it becomes stretched along the unstable direction, but still 
localized well enough to fit in the hole. It is therefore an approximate generalized 
eigenstate of the open map with vanishingly small eigenvalue. From |42| . reprinted 
with permission. 

U (the number of quantum states that fit in the hole) . Let A and 7 be the Lyapounov 
exponent and the escape rate of the system and let r e = i- In K be the (escape) 
Ehrenfest time. Remember that lZ~ m is the set of points that escape in exactly m 
steps. A minimal wavepacket supported in lZ~ m will therefore reach the hole in 
to steps. If m is smaller than r e , this wavepacket will not have spread too much and 
will still fit in the hole. It will then be almost completely annihilated and we will have 
/7 m |i/') ~ 0, i.e. it will be a generalized eigenvector of U with zero eigenvalue. This 
process is illustrated in Figure 5. 

What we mean by saying that \ip) is a generalized eigenvector is that and 
its images under U form a Jordan chain. In other words, the fractal Weyl law is 
related to the difference between the algebraic and the geometric multiplicity of the 
null eigenvalue. These ideas were further discussed in [74], where it was suggested 
that, instead ofjjeing diagonalizcd, the open map should be subject to a Schur 
decomposition, U — QTQ\ where Q is unitary and T is upper triangular. This 
has the advantage that, unlike the eigenvectors, the columns of Q are orthogonal. 

In essence, this identifies the states comprising the short-lived sector of the 
resonance spectrum. To estimate their fraction, we need only estimate the total area 
of all Hr rn with m up to r e . Since the fraction of points remaining after t steps 
is e~ 7 *, that area behaves like JI° e~ ll dt ~ 1- e ,Tc . Therefore, the short-lived 
sector amounts to a fraction of the spectrum that scales as N(l — if ~ 7 / A ). Since K is 
proportional to N, this implies that the total number of states in the long-lived sector, 
which is counted by the fractal Weyl law, should scale as N dl where d\ = 1 — 7/ A is 
the partial information dimension of the natural measure fi e . 

The above arguments provide some understanding of the physical origins of the 
FWL, but the exponent they predict is di(fi e ) instead of do(/j, e ). Rigorous upper 
bounds involve do- Still, the information dimension was also used in [71], where the 
kicked rotator is studied in a range of decay rates, and again in [72] for the cat map 
with several different holes. In both these works very good numerical agreement was 
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obtained. On the other hand, it was shown in |75j that the true exponent really is 
closer to 8q by explicit construction of a system where do and d\ are very different. In 
general, however, these two dimensions are quite close; large deviations can only be 
expected for systems where the local stretching factor is very far from being constant 
(which is the case in [75]). 

Another approach to the FWL was developed in [76] and [77]. This is based on 
semiclassical approximations for the long-lived states, and the basic idea is that these 
should be related to periodic points, which belong to the saddle. An approximate 
basis for the long-lived sector is constructed by building wavepackets concentrated 
on periodic orbits of period T, with T up to te = j\nN. These wavepackets are 
then propagated with the open quantum map, thereby acquiring some information 
about short-time dynamics and providing an improved basis, the so-called scar 
functions. Since the number of points with period less or equal to te grows like 
e (^-7)TE ^he quantity A — 7 is the Kolmogorov-Sinai entropy), the dimension of the 
basis is proportional to N dl . A matrix of this dimension is then introduced as an 
approximation of the quantum propagator in the long-lived sector. Its spectrum is in 
good agreement with exact results, validating the approximation. This is shown for 
the baker map in Figure 6. Notice the presence of a gap in the bulk of the spectrum. 

This semiclassical approximation was critically discussed in [7SJ , where it is argued 
that it reproduces the spectrum, but in order to correctly describe quantum resonance 
wave functions semiclassically it is necessary to take into account more information 
than only the chaotic saddle. The authors of [78] speculate that diffraction effects 
and/or classical information outside the saddle should be important and must be 
incorporated. This is yet to be systematically investigated. Another point raised in 
[75] is that perhaps in practice the fractal Weyl law only holds once h becomes much 
smaller than the size of the hole, otherwise it may be affected by diffraction. 

In |73) a rigorous theory was developed to the construction, based on an open 
quantum map of dimension A, of an auxiliary operator whose rank is ~ N d ° and whose 
spectrum reproduces the exact resonances. That operator is conceptually similar to 
the one based on short periodic orbits that was discussed above. 

Notice that the two approaches we just delineated are different and somehow 
complementary to each other. While [70] arrives at the FWL by estimating the size 
of the short-lived sector, associated with regions of fast escape, [7j3 [77] attempt to 
estimate the size of the long-lived sector, expected to be related to periodic orbits. In 
both cases some notion of Ehrenfest time plays an important role. 

The Weyl law for systems whose dynamics is a mixture of chaotic and regular 
regions was discussed in [74] . It was found that stability islands led to resonances with 
very small decay rate, approximately following a usual Weyl law, while states with 
intermediate decay rates obeyed a Weyl law with a fractional exponent; however, 
this exponent was not related to the dimension of any fractal set. The subject 
was also investigated in [75] for a system with sharply divided phase space. States 
associated with sticky motion at the border of the stability island satisfied Weyl laws 
with fractional exponents, related to the power law behaviour of the classical survival 
probability. 

5. Resonance gap 

The fractal Weyl law is a prediction about how the number of long-lived resonances 
scales with h or, for quantum maps, with the dimension N. A different type of 
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Figure 6. (color online) Exact spectrum of the open triadic baker map (circles) 
and the approximation to the long-lived sector obtained using scar functions 
(crosses). Four values of N are shown: 81, 177, 597, 1821. In each case, the 
approximation uses a basis of size proportional to A r_d / 2 where d is the dimension 
of Kq. Good agreement is obtained, validating the FWL scaling. The internal gap 
has not been explained. Figure reprinted with permission from |77| . Copyright 
(2012) by the American Physical Society. 

question is how these resonances are distributed in the complex plane. For maps, we 
denote them by z n = e * e "e~ r "/ 2 , so that \z n \ 2 = e~ r '\ In particular, it is natural to 
focus interest on the longest-lived ones because they usually leave clear signatures on 
scattering signals. Numerical experiments reveal that the distribution of decay rates 
r„ tends to have a maximum around the classical decay rate 7, with a long tail for 
T > 7 and a rather short one for T < 7. This has been observed for example in 
[7T1 l80l HU [82], and we show an example in Figure 7. 

In a sense, states with T < 7 have anomalously slow decay. They were termed 
'supersharp' resonances in [83 . A natural question is whether there is a limit to how 
sharp a resonance can be. Can the quantum decay rate be arbitrarily small? Evidence 
is on the contrary. It is generally expected that supersharp resonances should become 
less frequent in the semiclassical limit, i.e. that a true gap develops in the spectrum 
as N —> 00. 

Already in [84] (see also [85]) it was found that a lower bound could be proved 
for the decay rates of chaotic resonances in the limit H — > 0. Its formulation 
involves something called the topological pressure P(/3). This is related to the leading 
eigenvalue of a generalization of the Perron-Frobenius operator introduced by Ruellc. 
We do not discuss this so-called thermodynamic formalism here, and instead refer the 
reader to [86| HFJ. The topological pressure has some important properties, namely 
its value at 1 gives the decay rate, P(l) = —7, its derivative at this point gives the 
Lyapounov exponent, P'(l) = —A, and it has a zero at the partial dimension of the 
unstable manifold of the saddle, P(d(K-)) = 0. For the open triadic baker map it is 
given by P(/3) = In 2 — j3 In 3; in general it is a convex decreasing function. 

The available rigorous lower bound is that all r„ must be larger than — 2P(l/2) as 
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Figure 7. (color online) Distribution of decay rates, scaled by the classical one 7 C , 
of the open kicked rotator, for many different values of N. A limiting behaviour 
as N — ¥ oo and a sharp peak very close to 7 C are clear. Figure reprinted with 
permission from 71 . Copyright (2008) by the American Physical Society. 

N — > oo. This result has been revisited recently with more rigor in [57]- Interestingly, 
this bound is only effective if d(K^) < 1/2, because P(/3) is a decreasing function 
with a zero at d(K-), so P(l/2) > if d{KJ) > 1/2. This would just predict that all 
r„ are larger than a negative value, something that is always true. In other words, 
only for saddles that have low enough dimension, or are 'filamentary' enough, can a 
true lower bound really be established. 

It is known that 2P(/3) > P(2(3) and hence — 2P(l/2) is always smaller than the 
classical decay rate 7 = — P(l). Therefore, this lower bound still allows for quantum 
decay that is slower than classical. At present, it is not clear whether in the limit 
N — >• 00 one can still find resonances with decay rate arbitrarily close to — 2P(l/2) 
or if they tend to be larger than 7. In other words, if a larger gap can appear. An 
approach based on random matrix theory would predict the latter, as we will see in 
Section 8. 

Numerical experiments conducted for the baker map were somewhat inconclusive 
in this respect [25] . In this system, the smallest decay rate T seems to be influenced by 
the presence of discontinuity lines and seems to be different for even/odd eigenstates 
when the saddle intersects those lines. For the kicked rotator [55], on the other hand, 
it was clearly observed that e _r ° — e -7 approaches zero as a power law N~ s . In fact, 
the exponent S was found to be approximately given by d ^ ' — a, i.e. the exponent 
from the fractal Weyl law minus a constant. 

Also studied numerically in |83] was the distribution of supersharp resonances. 
It was found that they do not conform to the fractal Weyl law, i.e. their number 
does not increase in proportion to the bulk of the spectrum. Instead it seems that 
are approximately N a such states, where a is the same constant as in the previous 
paragraph. Whether this exponent is universal or system-specific and what classical 
information it contains are open problems. Notice that therefore the number of 
resonances with T < 7 grows with N, but at the same time the distance To — 7 
decreases with N. Similar results and conjectures have been surfaced in the study of 
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Figure 8. (color online) Left: two families of open baker maps that have different 
holes and different short-time dynamics, but the same decay rate and the same 
saddle. Right: for the 'intersection' family, the spectral radius decreases with k. 
Figure reprinted with permission from |90| . Copyright (2012) by the American 
Physical Society. 



the Laplace resonance spectrum on hyperbolic surfaces of infinite volume [HH US] • 

The smallest decay rate, unlike the Weyl law exponent, seems to be very sensitive 
to the short-time dynamics of the system, as observed in [3D]. In this paper two 
families of baker maps were defined, along with their quantizations, the shift and the 
intersection families, B^' and B^ for k > 1 (see left panel of Figure 8). Inside each 
family, all members have the same Lyapounov exponent and decay rate (and hence 
the same partial information dimension for their natural measures), but very different 
short time dynamics. In the intersection family, members have holes that grow with 
k, in such a way that only (2/3) 2k of all initial conditions remain after the first step. 
It was observed that the smallest decay rate grows with k, so that all resonances have 
fast decay for large k (see right panel of Figure 8). On the other hand, in the shift 
family the area of the hole is independent of k and so is Tq. 

6. Eigenvectors 

As discussed in Section 3, open quantum propagators are not unitary matrices and 
their left and right eigenvectors arc different, 

U\**} = z n \**), (^\U = z n {^\. (25) 

Just like for closed systems, one would expect on general grounds that in the 
semiclassical limit these wavefunctions should be related to classical structures. A 
natural question is the analogue of the quantum ergodicity problem: what are the 
possible semiclassical limits of Husimi functions of resonance wavefunctions? 

The issue of quantum ergodicity for open maps is actually much more difficult 
than for closed ones, because there are many conditionally invariant measures for any 
given decay rate, and it is not clear which ones would be more natural as semiclassical 
limits of Husimi functions. 

A few general results were obtained in 91j. For example, Husimi functions of 
right eigenstates, H^R(q,p) — \(q,p\^>^;)\ 2 , become concentrated, in the semiclassical 
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Figure 9. Two Husimi function of right eigenstates of the open triadic baker 
map, for N = 1500. Intensity grows from light to dark. Eigenvalue modulus is 
0.84 for the left one, which is approximately uniformly distributed along A*_ , and 
0.46 for the right one, which shows some localization on the hole and 7£ _1 . From 
25 , reprinted with permission. 





Figure 10. (color online) The weight of Husimi functions of right eigenstates of 
the open triadic baker map on the regions 1Z~ m , for N = 729. The solid line is 
the semiclassical approximation d 26 1) . Figure reprinted with permission from |91| . 
Copyright (2006) by the American Physical Society. 



limit, on the backward trapped set K~. This is true in the following sense: suppose 
a sequence of resonances with increasing N, such that the corresponding eigenvalues 
converge to some value z which is different from zero. Then the sequence of (right) 
Husimi functions will converge to zero at points which do not belong to K _ . We show 
two illustrative examples in Figure 9. Analogously, left eigenstates concentrate on K + 
as N — > oo. The results in [21] do not indicate how fast the value of 'H^R{q,p) goes to 
zero if (q, p) ^ K_ , but in [75] it was shown that this happens (at least for the baker 
map) exponentially fast, i.e. like e~ cN for some constant c. 

It is also possible to obtain some semiclassical information about how these 
Husimi functions are distributed on their support. First, remember that the sets lZ~ m 
intersect the backward trapped set. Let these regions be semiclassically quantized by 
phase space projectors ir m (as discussed in |M]). The open propagator thus satisfies 
U'U = 1 — 7To . It was shown in [HI] that 

(*%\n m \*%)*\z n \ 2m (l-\z n \ 2 ). (26) 

The left hand side of the above equation measures the weight of 7-L^R{q,p) on the 
region lZ~ m , and the right hand side relates this to the decay rate. This relation is 
tested in Figure 10 for the baker map. 
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It was first noticed that the functions 7i^,R(q,p) had fractal signatures of the set 
K- in [55]. In this paper it was suggested that, for a sequence of states whose decay 
rate converges to the classical one, the (right) Husimi functions should converge to 
the classical equilibrium measure. This would in fact be consistent with the result 
(|26|) above, because since /x e (/i) = 1 — e~ 7 and since the area decreases by e -7 at each 
step, we have 

fi e (TZ- m ) = e- m7 (l - e~ 7 ). (27) 

However, this convergence has not been proven and numerical evidence is inconclusive. 
Equation (|21)|) also predicts that states with r„ > 7 must show some concentration on 
the hole and its first pre-images and, conversely, that states with T n < 7 (supersharp) 
must avoid these sets and localize closer to the invariant set. 

This problem was extensively discussed for the baker map in [75] . In particular, 
this system admits a non-canonical quantization, known as the Walsh quantization, 
in which it is exactly solvable (unfortunately it is not clear how general are the results 
proven in this special setting) . For example, right Husimi functions having self-similar 
properties can be explicitly constructed for this system. This was discussed again in a 
slightly different setting in |93j where, among other things, a kind of quantum unique 
ergodicity was proven at the edges of the spectrum. 

Related to the problem of quantum ergodicity is the phenomenon of quantum 
scarring. This is a generic name given to concentration of quantum states in the 
vicinity of periodic orbits. For open systems the question of scarring carries an extra 
interest, because of the interplay with the decay rate. States with small decay rate 
must survive in the system for a long time, and it is natural to expect them to be 
scarred along periodic orbits. States with T > 7 can show scarring, and this has 
been observed. But, as discussed in relation with (|26j) . scarring is more likely for 
T < 7 (this was also suggested in [81]). Moreover, how the effect behaves in the 
semiclassical limit is rather unclear. This is connected to the subject of the previous 
Section, because a gap may develop in the spectrum such that no states with T < 7 
survive the semiclassical limit. 

It was observed in [82] that resonances of the kicked rotator with small decay rates 
were scarred. The subject was taken up again in [72], where the open cat map was 
studied. A set of quantum states specially suited for studying scarring, the so-called 
scar functions already available for closed maps [MJ [SS], was adapted to the open 
setting in [761 [77] , where it was observed that the baker map with low N had some 
resonance wave functions that could be obtained as a superposition of only two of 
these states. However, it was also suggested in [78] that classical information outside 
the saddle, such as diffraction effects, must somehow be incorporated if a semiclassical 
approach is to accurately reproduce resonance wave functions. 

In this context, a new phase space representation of resonance states was 
introduced in |96] , which is a generalization of the Husimi function. It is given by 

Since {q,p\^^) is localized on K- and (^^\q,p) is localized on K + , the quantity 
h n (q,p) must be localized on their intersection, the saddle, and may be useful for 
revealing scarring effects. This was verified for the baker map. It was also noticed 
that, for finite N, the states with large decay rates can have significant values of h(q,p) 
outside of that fractal set. 
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Figure 11. The scarring estimator r\ for the open triadic baker map at N = 243. 
The circles correspond to individual states, rj(hj), ordered with increasing decay 

rate. The line corresponds to the cumulative value r\ ^^Zi=i hi). Courtesy of L. 
Ermann. 



where 



A quantity derived from h(q,p) was introduced in [90] as a quantitative measure 
of localization: 

(h) f \h\i/\h\ 2 \ 2 

\h(q,p)\ s dqdp) (30) 

and p represents the density matrix of any coherent state. As defined, r\ is equal to 1 
for a coherent state (perfect localization) and N/2 for a uniformly distributed state in 
a A^-dimensional Hilbert space. We see in Figure 11 that the value of r\ behaves rather 
erratically as a function of the decay rate for individual states. However, a cumulative 

value r\ (^J2i=i h-ip where the sum runs over states with increasing decay rate, shows 
an almost monotonic behaviour, with smaller decay rates corresponding to stronger 
localization. 



7. Non-ideal escape 

A generalization of open systems which has attracted attention recently is the 
possibility to define holes that are only partially transparent: when a point hits the 
hole, it has a finite probability of escaping, but may also continue being propagated 
inside. This may be called non-ideal or refractive escape, because it is supposed to 
mimic the ray-splitting that takes place at the interface between two media with 
different indices of refraction. It can also be used to model quantum dots coupled 
with tunnel barriers or any kind of losses. 

Let us suppose a dielectric sample of refractive index n, outside of which the 
refractive index is 1. When the wavelength is much smaller than the average radius of 
the sample, we may consider propagation of rays. These follow straight lines, except 
when they meet the boundary. According to Snell's law, they then experience specular 
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reflection. This is what is known as a dynamical billiard. It is usual to define boundary 
coordinates as follows. Choose an arbitrary point r at the boundary and let q be the 
distance along the boundary between a collision and r . Clearly, < q < L, where 
L is the perimeter of the sample. The other coordinate, p, is the sine of the angle 
between the incoming ray and the normal to the boundary at the collision point. It 
is possible to show that in this coordinate system the dynamics is area-preserving. 

The propagation of rays inside the dielectric sample is thus reduced to a discrete 
map in the variables (q,p), of the kind we have been considering. Now we must 
incorporate the fact that, at each collision, some amount of light is refracted out of 
the sample. The intensity of the reflected beam is given by Fresnel's law (for transverse 
magnetic polarization, for instance): 



and depends on the angle of incidence. Refraction does not occur if \p\ > p c = 1/rt, 
which is called the critical angle. 

Refraction can be modeled in a chaotic map by assigning an escape probability 
to each point in phase space [97]. An attempt to be realistic would try to simulate 
Fresnel's law. More simply, one can define a finite hole corresponding to \p\ < 1/n and 
associate to it a constant escape probability. Now it is no longer possible to simply 
follow initial conditions with time, because one must keep track of probabilities. This 
can be solved by considering only the evolution of probability densities (this is also 
discussed in [22]). 

We can start with the uniform density in phase space, for example, and let it 
evolve under the partially open map. Equivalently, we can associate an initial intensity 
with each point, and this intensity is reduced every time the point falls in the (partially 
transparent) hole. In particular, the set of points whose intensity remains forever equal 
to its initial value is nothing but the forward trapped set that would correspond to a 
totally transparent hole. However, other regions may also have high intensity and be 
important for long-time properties. 

Strictly speaking, the system will have a (conditionally invariant) equilibrium 
measure [i e whose support (the 'forward trapped set') is the entire phase space, since 
the intensity of a point will never be exactly zero. In other words, pL e will have a 
trivial Minkowski dimension do(/i e ) = 2. However, the intensity after a long time will 
be a wildly fluctuating quantity and usually there will exist non-trivial dimensions 
d/s(fi e ) < 2 (naturally, the same is true for backward time evolution). Further 
discussion of these points can be found in [98] [99] [22], along with pictures of such 
intensities for the stadium, annular and cardioid billiards, respectively. We show an 
example in Figure 12. 

Quantization of partially open maps poses no difficulty. We continue assuming 
that the hole is a strip (or a union of strips) parallel to one of the axes. Instead 
of using a projector n, one needs only multiply the closed-system propagator U by 
a diagonal matrix whose elements are equal to ^/R where the reflection probability 
is R (in principle, the phases of such elements are arbitrary; obviously, we assume 



The first question one would ask in this setting is about the analogue of the Weyl 
law. Do we see a fractional exponent in the way the number of resonances scales 
with Nl It was initially suggested in [98] that the Weyl law could be sensitive to the 
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Figure 12. (color online) The conditionally invariant natural measure of the 
dielectric annular billiard, with refractive index equal to 4, depicted in the 
boundary coordinates described in the text (vertical axes is p, horizontal one 
is q). Intensity grows from light to dark. Figure reprinted with permission from 
1 991 . Copyright (2009) by the American Physical Society. 



different fractal dimensions present in the classical equilibrium measure. By counting 
resonances in the dielectric stadium billiard, different exponents in the Weyl law were 
found depending on the decay rate; moreover, these exponents were in the same range 
as the numerically determined dp(fi e ). 

However, it was later proved jlOOl 1101] that in the semiclassical limit the 
resonances tend to cluster at a particular value of the decay rate. More specifically, 
let R(q,p) denote the reflection probability at point (?,£>)■ In the long run, a typical 
(ergodic) trajectory will sample all of phase space. At each point, its intensity is 
multiplied by R. This process is additive in the logarithm, resulting in the average 
decay e~^ r \ where 



It was shown in jlOOl 1101] that the number of resonances with decay rate inside a 
fixed interval containing (r) scales with N 1 in the semiclassical limit. In other words, 
most resonances have this decay rate. This result is illustrated in Figure 13. Since the 
resonances counted in [95] were in a range that contained (r) , the fractional exponents 
observed are probably due to insufficient statistics or the semiclassical regime not been 
reached yet. 

Another result presented in [100] 1101] is that, for strongly chaotic systems, the 
width of the distribution of decay rates around (r) decays like (IniV) -1 . This was 
proven for a toy model and numerically observed (albeit not quite clearly) for the baker 
map and a perturbed cat map. It is natural to expect that this spectral width should 
depend on some type of variance of the function R(q,p). In particular, the width is 
zero if R is constant all over phase space. This point has not been investigated. 

The distribution of resonances for partially open maps had already been 
considered numerically in [57] using the kicked rotator as a model. This paper also 
suggested a random matrix theory approach (see next Section), which was later tested 




(32) 
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Figure 13. (color online) Cumulative distribution function of decay rates for 
partially open baker maps (curves a\ and ci2 correspond to two different choices 
for hole position and reflectivity). Vertical lines correspond to (r), and clearly 
the distribution has an inflection point there. Notice how little change there is 
with N. Figure reprinted with permission from 11001 , Copyright (2008) by the 
American Physical Society. 



in [102] for the stadium billiard. Good agreement was obtained, as long as the mean 
decay rate and the effective number of modes were suitably adjusted. 

The phase space morphology of resonance states in the case of non-ideal escape 
has not been investigated for maps, as far as this author is aware (see however 

doa una una). 

8. Random matrix theory 

An approach which has been very successful in all areas of quantum chaos is the theory 
of random matrices [106j . In essence, the specific relevant operator (Hamiltonian, 
S matrix, propagator) is abandoned and replaced by a random matrix, i.e. only 
statistical properties are studied with respect to a certain ensemble of matrices. How 
this ensemble is chosen is dictated by the situation at hand. Approaches based on 
random matrices exist for several areas of physics, including the quantum mechanics 
of closed chaotic systems |107] and chaotic transport [28]. Many examples can be 
found in a recent Handbook [108] , 

To replace the Hamiltonian of a closed system, for example, it is important that 
hcrmitian matrices be used. The simplest possibility is then to replace the matrix 
elements with identically distributed Gaussian random variables, while respecting 
the hermiticity constraint. This gives rise to the Gaussian Unitary Ensemble (it is 
called Unitary because it is invariant under unitary transformations). If there is 
time-reversal symmetry, the Hamiltonian can be made real without diagonalization. 
What is needed then is an ensemble of real symmetric matrices, and this leads to the 
Gaussian Orthogonal Ensemble. It is common practice to associate an index (3 = 2 
with complex matrices and j3 = 1 with real matrices (this is sometimes called the 
Dyson parameter). Gaussian ensembles were recently reviewed in |109j . 

8.1. Effective Hamiltonian approach 

It is possible to model scattering by writing the S matrix as S(E) = 1 — iV^ E _j { ^ V, 
where V is a N x K coupling matrix between the system and the outside and H e tf 
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is an N x N effective non-hermitian Hamiltonian. The poies of S are then the 
eigenvalues of H e fj. This is sometimes called the Heidelberg approach. Assuming 
energy-independence of the coupling elements and neglecting direct processes, this 
latter matrix can be taken as H e ff = H — \VV* , and H considered as being 
uniformly distributed in one of the Gaussian ensembles. The reader can consult 
PUJ E21 ESI ES [24] for further information. 

Different regimes can be studied depending on coupling strength and the interplay 
between K and N. The limit N — > oo is always intended. The number of channels K 
may remain fixed, in which case the hole is not classical in size. For small coupling, 
all eigenvalues acquire small imaginary parts (smaller than the mean level spacing). 
For large coupling, the phenomenon of resonance trapping occurs |114) : only a few 
resonances acquire large imaginary parts, while most of them approach the real axis. 

The case of classical opening, closer in spirit to the subject of this review, 
corresponds to keeping the ratio K/N fixed as N — > oo. This corresponds to the 
size of the hole, K/N = h. As we have seen, it makes sense to consider this as a small 
variable. It can be shown |115[ 1116] that in this case the density of resonances exhibits 
a gap: it vanishes for decay rates smaller than h, which can be identified with the 
classical decay rate. Moreover, this density decays like 1/r 2 for V > h. 

8.2. Approach via truncations 

The propagator of a closed quantum map is a A-dimensional unitary matrix, UU^ = 1. 
If there is time-reversal symmetry, this matrix is also symmetric. These are the only 
constraints that must be imposed. The ensemble of random unitary matrices is nothing 
but the unitary group equipped with its unique invariant (probability) measure, the 
Haar measure. In this context, this ensemble is called the Circular Unitary Ensemble, 
the CUE. If symmetry is imposed, U T — U, we obtain the Circular Orthogonal 
Ensemble, the COE. 

The propagator of an open quantum map is obtained by multiplying the closed 
propagator by a projector. This is sometimes called a truncation: some columns of 
U are set equal to zero to produce U. Therefore, it is natural to take U as being 
uniformly distributed in the CUE and introduce an ensemble of truncated unitary 
matrices (TCUE). This was done in [TTO] . 

We therefore have iV-dimensional unitary matrices truncated by multiplying from 
the right by a projector. Let the kernel of the projector have dimension K, as before, 
and let M — N — K be the number of eigenvalues that are not identically zero. Just 
like for the effective hamiltonian approach, in the 'semiclassical' limit N — > oo there 
are two choices for the behaviour of the truncation: the kernel dimension K can remain 
fixed in order to model 'quantum' holes, or we can have M = jj,N with fixed \x in order 
to model 'classical' holes. 

Since any element of the TCUE is subunitary, the spectrum always lies inside the 
unit disk in the complex plane. The joint probability density of the eigenvalues can be 
exactly obtained. For large N, the original ensemble of U is not relevant, i.e. results 
are the same for both unitary and orthogonal classes. Apparently, the presence of 
the hole is enough to effectively break symmetry. The probability distribution of the 
modulus squared of the eigenvalues can also be found. Using the quantity /i mentioned 
above, the distribution of p = \z\ 2 = e~ r is given by 

l-H 1 
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if p < p and vanishes otherwise. 

Since p = 1 — ft, is the area of the complement of the hole, it is natural to interpret 
it in terms of a fictitious decay rate as p — e~ 7 . Therefore, this approach also predicts 
a spectral gap in the limit N — > oo, so that states with T < 7 become increasingly rare. 
Note however that the interpretation of T as a decay rate is not so straightforward for 
discrete time systems (maps), e.g. 

e -rt/h dt= n wh . le y e -m = — L (34) 
r ^ 1 - e~ r v ' 

n=0 

If we were to interpret 1 — e~ r as a decay rate, then according to ()33|) its density 
decays quadratically, in agreement with the effective Hamiltonian approach. 

In Figure 14 the prediction (|33p is checked for the open kicked rotator [not 
available in the Arxiv version]; good agreement is found, provided the dimension 
N is renormalized according to the fractal Weyl law. 

Notice that this approach is not able to give any insight into the fractal Weyl law 
itself, since this is obviously a system-specific feature, i.e. it depends on details of the 
dynamics. One way to see this is that although it is possible to build the escape rate 
into random matrix theory, simply using the size of the hole, the same is not true for 
the Lyapounov exponent (in a sense, this quantity is 'infinite' in this theory). 

Even though they are suppressed in the N — ¥ 00 limit, we can ask about statistics 
of resonances inside the gap, the supersharp ones. For example, we can consider states 
with p = p, + e/77, with i] = \J N/2p(l — p). These will approach the maximum 
value from above in the semiclassical limit. It turns out that they have a finite 
density function proportional to erfc(e), where erfc is the complementary error function 
|117j . Another interesting question is the distribution of the smallest decay rate Tq. 
Equivalently, we can consider the largest eigenvalue modulus squared, po = e~ r °. It 
was shown in |118l 185] that the modified variable rj 2 (p — p) 2 — In N + In In N satisfies 
a Gumbel distribution as N — > 00. 



L 



8.3. Eigenvectors 

A different line of investigations would be to consider statistics of the eigenvectors. 
Usually, the quantity 

(\fr L \\S) L )(\fr R \\lr R ) 

tv- \ I nl n/\ ?i I n/ / o c: A 

- " mm? ' [ J 

is considered (see, for example, (1191 1120] ). This is called eigenvector correlator 
or Petermann factor, and is related to the line width of the lasing mode in open 
resonators. Results of this type were reviewed in |112) . This can be generalized to an 
off-diagonal version U nm = (^n\^m)(^n\^m)> which was shown in |121] to be related 
to the statistics of resonance width shifts under external perturbations. 
Another type of quantity, 

where ^ is the ith component of a right eigenvector, was considered in [1221 1123] . 

Research along this line has always been done within the effective Hamiltonian 
approach, in the case when the number of states in the hole is much smaller than total 
number of states, K <gC N. To the knowledge of this author, there are no works about 
statistics of eigenvectors for classical holes or truncated unitary matrices. 
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8.4- Non-ideal escape 

The case of non-ideal escape can also be treated in two ways, using effective 
Hamiltonians or opening propagators. Within the first approach the finite 
transparencies of the decay channels are incorporated into the coupling matrix V. 
When all channels are equivalent the degree of resonance overlapping is controlled by 
the parameter KT |124j where T is the transparency and K is the number of channels. 

The propagator of partially open systems should be modeled as U = U VG 7 
where G is a fixed diagonal matrix with nonnegative entries < g\ < gi < • • • < 1 
representing reflection probabilities and U is uniformly distributed in a circular 
ensemble. This was first considered in [97] . with U in the COE, but only numerical 
results (for distribution of decay rates and eigenvectors) were presented. Later, the 
exact density of states with U in the CUE was derived in [125] . Its asymptotic 
behaviour as N — > oo, obtained in |1261 1127] . is as follows: if x — \z\ denotes eigenvalue 
modulus, its cumulative distribution y(x) is implicitly defined by 



In particular, in the semiclassical limit the eigenvalues of U lie inside an annulus in 
the complex plane: their modulus squared must be larger than the arithmetic mean 
of the eigenvalues of G and smaller than their harmonic mean. 

9. Some open problems 

The topic of resonances is of wide interest, both theoretically and experimentally. 
There is a variety of settings, motivations, applications, approaches, etc. that make 
it impossible to provide a truly comprehensive review. We have focused on quantum 
maps, with chaotic classical dynamics and classical openings. We hope to have touched 
upon some of the more interesting points that have recently attracted attention and 
to have provided some references for the interested reader. 

We conclude this review with a few open problems. These are more general 
indications of interest than specific questions. It is likely that as these matters are 
investigated further other interesting problems will come to the front. 

For systems with ballistic escape, we would like to mention the following points: 

• Exponent in the fractal Weyl law. A rigorous proof of this, aside from the 
non-generic Walsh quantization of the baker map, is still lacking. Heuristic 
understanding has been achieved from different points of view, however, and 
such proof is likely to be extremely technical. 

• Prefactor in the fractal Weyl law. Does the profile function of the resonance 
spectrum C (r) have anything to do with classical dynamics? The random matrix 
theory prediction works well for the kicked rotator, but not for the baker map. 
Can this be understood? 

The next points apply to ballistic escape, but also make sense for systems with 
non-ideal escape: 




(37) 



where \& is 




(38) 
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• Spectral gap. Under what conditions is it true that as N — > oo there are no 
resonances in a certain region? What exactly is the Weyl law for supersharp 
resonances? 

• Quantum ergodicity. What are the possible semiclassical limits for Husimi 
functions with a given decay rate? Do they have self-similar properties? Is there 
a preferred classical eigenmeasure? 

• Scarring. How does the scarring by periodic orbits behaves in the semiclassical 
limit? Is the amount of scarring related to the decay rate? Are supersharp 
resonances special in this respect? 

• Semiclassical approach. It seems that in order to reproduce resonance wave 
functions it is necessary to include more information than just periodic orbits 
(maybe diffraction). It is not clear how to proceed in this direction. 

• Statistics of eigenvectors. Virtually nothing is known in the case of classical 
openings, even within random matrix theory. 
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